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Abstract 

' Method of multiple internal reflections (method MIR) in description of a-decay of nucleus in the spher- 

ically symmetric approximation is presented in paper. In approach MIR the formalism of calculation of 
amplitudes of wave function, described moving of a-particle from the internal region outside with its tun- 
neling through a realistic radial a-nucleus barrier of arbitrary shape, has been constructed at first time. 
The method MIR gives convergent values for the amplitudes of transmission and reflection, coefficients of 
^ . penetrability and reflection, obtained in description of leaving of the Q-particle relatively a potential with 

^ ' barrier in form of a number of rectangular steps (multi-steps potential), with tending this potential to the 

, realistic a-nucleus one, and limit values of the amplitudes and coefficients are considered as exact concerning 

the realistic a-nucleus potential (without application of WKB approximation) . In approach MIR a depen- 
dence of the penetrability coefficient and half-live on the location of the starting point, from where a-particle 
begins to move outside, is opened. From here, we define a radial coordinate of the most probable formation 
1^-^ ' of the Q-particle inside nucleus before its a-decay as such starting point, at moving of the Q-particle from 

\^ I which outside the half-live, calculated by approach MIR, is maximally closed to its experimental value. For 

isotopes Po with different mass number A we have obtained essentially closer values for half-live, calculated 
by such approach, to their experimental values in a comparison with half-lives obtained by approach WKB. 
, This result is stable for all studied nuclei. 

O ; PACS numbers: 03.65.Xp, 23.60.+e, 27.80. -hw 

Keywords: tunneling, multiple internal reflections, wave packet, alpha-decay, coefficients of penetrability 
and reflection, half-live, times 
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1 Introduction 

An approach for description of one-dimensional motion of a non-relativistic particle above a barrier on the 
basis of consideration of multiple internal reflections of stationary waves relatively boundaries of this barrier, 
described a motion of this particle inside the barrier region, was considered in a number of papers and is known 
[TJ 1^ [3]. In such approach it needs in expression of wave function (WF) in the barrier region to separate 
components, having fluxes directed to opposite sides. For a problem of the motion of the particle above a 
rectangular barrier these components are found naturally — as plane waves e^^^^ where k is wave vector. To 
obtain solutions for description of tunneling of this particle under such barrier it turns out more difficult, because 
in consideration of tunneling as stationary process the fluxes defined on the basis of decreasing and increasing 
parts of the stationary WF in dependence on x inside the barrier region (which are analytical continuations of 
the corresponding waves for case of above-barrier propagation) equal to zero and it is not corrected to consider 
such "waves" as propagated (or tunneled) under the barrier from physical point of view. If to define the waves 
under the barrier by another way (for example, at first to require existence of non-zero fluxes in the barrier 
region in each step and then to find solutions of the waves), then we cannot obtain continuity in transition of 
expressions for the waves for the above-barrier and sub-barrier cases, and total amplitudes will be different on 
the corresponding amplitudes, determined by standard approach of quantum mechanics as, for example, in [3] 
(see p. 102-106; for simplicity we shall further name such approach for determination of the amplitudes of the 
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stationary wave function as standard (direct) method of quantum mechanics) . However, the flux, defined on the 
basis of the total stationary WF, is non-zero, that points out that tunneling of the particle under the barrier 
exists. So, there was a difficulty in correct description of tunneling of the particle under the barrier on the basis 
of multiple internal reflections of the stationary waves. 

However, this question was successfully resolved at first time in papers [SJ [S] (see also [7], part 1), in 
result of introduction of wave packets (WPs). In such approach, tunneling of the particle under the barrier is 
considered on the basis of multiple internal reflections of the wave packets relatively boundaries of the barrier 
(this approach was named as method of multiple internal reflections or method MIR). In such approach it 
succeeded in connecting: 

• continuous transition of solutions for separate WPs concerning each reflection from multiple ones (i. e. 
in each step, according to formalism of method MIR), total WPs (and corresponding stationary WFs) in 
transition from the above-barrier motion to the under-barrier tunneling and back; 

• coincidence of the transmitted and reflected total amplitudes of stationary part of WF (in each spacial 
region of the potential) in approach MIR with the corresponding amplitudes obtained by standard method 
of quantum mechanics; 

• all fluxes, constructed on the basis of WPs in each step, are non-zero, that allows to consider propagation 
of multiple packets under the barrier (i. e. their "tunneling"). 

As a further development of a non-stationary formalism for description of tunneling processes in direction of 
papers [H [9], this method has own description (interpretation) of non-stationary tunneling of the particle, 
allowing to study this process at interesting time moment or space point. In calculation of phase times this 
method has shown as enough simple and convenient p^. The formalism of MIR was developed also for 
description of scattering of the particle on nucleus and a-decay in the spherically symmetric approximation 
with the simplest radial barriers in [5l [5] (see also [7], part 2), and for tunneling of photons (for example, 
see [SmU]). 

However, in papers [1] [2l [3l [6l [10] correctness of applicability of the multiple reflections was proved for 
description of the motion of the particle above the barrier on the basis of waves, and also for description of 
tunneling of the particle under the barrier on the basis of wave packets, if this barrier has a shape of one 
rectangular step. But, addition of the second step to such barrier complicates essentially consideration of 
the multiple reflections of packets (or waves) relatively all boundaries and search of exact solutions for the 
amplitudes. It becomes unclear how to separate the needed internal reflections from all their variety inside 
interesting spacial part of the barrier for determination of the amplitudes here (for example, if the barrier 
consists of arbitrary number of the rectangular steps with different heights). So, a necessity in development of 
clear algorithm of calculation of all amplitudes in problems with the barriers of complicated shape has been 
appeared. 

In [5] exact solution of the amplitudes of the stationary wave function for the barrier, located in the whole 
axis and composed of two rectangular steps of arbitrary shapes, was found. However, there was a question in 
determination of wave packets under the barrier on the first step (this was resolved in 7j, part 1) and it was 
uncleared how to generalize such a formalism after complication of the barrier shape. Some late, in [llj approach 
of multiple internal reflections of the waves was considered for tunneling of the particle through a number of equal 
rectangular steps located in whole axis and separated on equal distances one from another. However, solutions 
of the amplitudes, the coefficients of penetrability and reflections were presented for two such steps only, in a 
approximation when they are separated on enough large distance. Besides, the found solutions in approach of 
multiple internal reflections are constructed on the basis of the amplitudes of total wave function, obtained early 
by standard direct method (see Appendix A, (7), (18), (19) in [H]). Therefore, a question about convergence of 
sums of the amplitudes of the multiple waves by approach MIR with the corresponding amplitudes of total wave 
function by standard direct method inside each spacial region (and, especially, inside the barrier region) remains 
open. So, we come to a serious problem of realization of the approach of multiple reflections in real quantum 
systems with barriers (solution of which can call in question effectiveness of applicability of the approach MIR, 
at all). 

After complication of the barrier shape the second question has been appeared: Whether is interference be- 
tween separate waves (packets) appeared, which are formed after separate reflections relatively different boundaries 
and reduces (or increases) the total wave? Whether does this come to principally different results of the approach 
of multiple internal reflections and direct method of quantum mechanics ? Note that such interference cannot be 
appeared in the problem of tunneling through one rectangular barrier (and motion above it) and, therefore, in 
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previous papers it was not visible. This problem has already been calling in question a correctness of approach 
MIR in application to the barriers with the complicated shapes. 

The given paper is directed on study of such questions and gives answer on them. In Sec. 2 the formalism 
of the method MIR in description of tunneling of the particle under one-dimensional rectangular barrier and its 
motion above it in whole axis (based on [Sl[5], [7] part 1) is presented. Here, we put proof of the method, analyze 
its peculiarities, write solutions for amplitudes (without phase times). Further, we present exact solutions of 
the amplitudes, the coefficients of penetrability and reflection by the method MIR in description of tunneling 
of the particle through arbitrary number of potential rectangular steps with arbitrary heights and widths in 
full consideration of the multiple internal reflections relatively all boundaries (at first time, a main idea was 
introduced in [7], see p. 127-129). According to analysis, the total amplitudes of transmission and reflection 
inside each region and concerning the total barrier coincide with the corresponding amplitudes obtained by 
the standard method. Such formalism transforms the method of MIR, early used mainly in problems with one- 
dimensional rectangular barriers, into a power tool of calculation of exact solutions of WF, the coefficients of 
penetrability and reflection relatively one- dimensional potentials of arbitrary shape, used in description of real 
quantum systems. 

In Sec. 3 we present the formalism of MIR in study of tunneling processes in problems of scattering of a par- 
ticle on nucleus and a-decay of the nucleus in the spherically symmetric approximation. At first, both processes 
are studied concerning a radial barrier of the simplest shape - as rectangular step. In the scattering problem 
the packets described evolution of this process (n-multiple WP formed in result of successive n transmissions 
and reflections relatively the boundaries of the barrier; total WPs inside each region; total WP transmitted 
through the barrier and total reflected WP from it), phase time of tunneling of the particle through the barrier 
and phase time of its reflection from the barrier are found. Expression for S-matrix is presented as a sum of 
two components, which are defined by the amplitudes of the stationary parts of the transmitted packet through 
the barrier and the refiected packet from it and correspond to resonant and potential scattering of the particle 
on the nucleus. 

Using idea in [12] , a main attention in this paper we concentrate to development of the formalism of MIR in 
description of the a-decay. At first, the a-decay with a radial rectangular barrier is considered. Here, the packet 
transmitted through the barrier which describes propagation of the particle from the internal region outside 
after its tunneling through the barrier, and the packet reflected from the barrier which describes delay of the 
particle inside the internal region in result of its reflection from the barrier are found; non-stationary analysis 
is fulfilled where duration of a-decay of nucleus is defined and calculated (see [5, part 2). Further, in Sec. 5 
we at first time present the formalism of MIR in description of the a-decay with a realistic a-nucleus radial 
barrier, when this barrier is approximated by arbitrary number of rectangular potential steps. In such approach, 
we demonstrate convergence of found solutions for the amplitudes of transmission and refiection for leaving of 
the a-particle relatively such multi-steps potential with barrier, tending such potential to the initial realistic 
a-nucleus one. Limit values of the amplitudes and coefficients are considered as exact concerning the realistic 
a-nucleus potential (without application of approximation WKB), on the basis of which we define half- live. In 
approach MIR a dependence of the penetrability coefficient and half-live on location of starting point i?start is 
observed, from where the a-particle begins to move outside. On this basis, we define a radial coordinate of the 
most probable formation of the a-particlc inside nucleus before its a-decay as such starting point, at leaving of 
the a-particle from which the half-live is maximally closed to its experimental value. It turns out, that half-lives 
calculated by such approach for a number of nuclei with the same Z = 84 and different mass number A, are 
closer essentially to their experimental values in a comparison with half-lives found in 13J by approach of WKB, 
that demonstrates effectiveness of the method MIR in description of the a-decay with the realistic a-nucleus 
barriers and in calculation of half-lives. 

2 Tunneling of the particle through barrier in one-dimensional prob- 
lem 

At first, we shall describe a general formalism of multiple internal reflections, using a simplest problem of 
tunneling of the particle through one-dimensional rectangular barrier in whole axis. 

2.1 Tunneling of wave through one rectangular barrier in whole axis 

Let's consider a problem of tunneling of a particle in a positive x-direction through an one-dimensional rectan- 
gular potential barrier (see Fig. [T|). Let's label a region I for a; < 0, a region II for < x < a and a region HI 
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for X > a, accordingly. Let's study an evolution of its tunneling through the barrier. 
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Figure 1: Tunneling of the particle through one-dimensional rectangular barrier 

In a beginning we consider a standard approach to a solution of this problem 0J[T3]. Let's consider a case 
when levels of energy lay under than a height of the barrier: E < Vi. The tunneling evolution of the particle 
can be described using a non-stationary consideration of a propagating WP 



where the stationary WF has a form: 

ip{x) = 



+00 

J g{E - E)ip{k, x)e-'^''^dE, 



(1) 
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for a: < 0; 
for < a; < a; 
for X > a: 



(2) 



and k = j-V2mE, ^ = j-y^2m{Vi - E), E and m are the total energy and mass of the particle, accordingly. 
The weight amplitude g{E — E) can be written in a form of gaussian [8j and satisfies to a requirement of the 
normalization J \g{E — E)\'^dE = 1, value E is an average energy of the particle. One can calculate coefficients 
At, Ar, a and f] analytically, using a requirements of a continuity of WF ip{x) and its derivative on each 
boundary of the barrier. 

Substituting in Eq. |T]) instead of ip{k,x) the incident tpincik, x), transmitted Lptr{k,x) or reflected part of 
WF ipref{k,x), defined by Eq. we receive the incident, transmitted or reflected WP, accordingly. 

We assume, that a time, for which the WP tunnels through the barrier, is enough small. So, the time 
necessary for a tunneling of an a-particle through a barrier of decay in a-decay of a nucleus, is about 10^^^ 
seconds. We consider, that one can neglect a spreading of the WP for this time. And a breadth of the 
WP appears essentially more narrow on a comparison with a barrier breadth [51 [S]. Considering only sub- 
barrier processes, we exclude a component of waves for above-barrier energies, having included the additional 
transformation 

g{E~E)~.giE-E)9iVi-E), (3) 
where ^-function satisfies to the requirement 

^ \ 1, for 77 > 0. 

The method of multiple internal reflections considers the propagation process of the WP describing a motion 
of the particle, sequentially on steps of its penetration in relation to each boundary of the barrier [T][2l[3]. Using 
this method, we find expressions for the transmitted and reflected WP in relation to the barrier. 

At the first step we consider the WP in the region I, which is incident upon the first (initial) boundary of 
the barrier. Let's assume, that this package transforms into the WP, transmitted through this boundary and 
tunneling further in the region II, and into the WP, reflected from the boundary and propagating back in the 
region I. Thus we consider, that the WP, tunneling in the region II, is not reached the second (final) boundary of 
the barrier because of a terminating velocity of its propagation, and consequently at this step we consider only 
two regions I and II. Because of physical reasons to construct an expression for this packet, we consider, that 
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its amplitude should decrease in a positive a;-direction. We use only one item /3exp(— ^x) in Eq. throwing 
the second increasing item aexp(^a;) (in an opposite case we break a requirement of a finiteness of the WF for 
an indefinitely wide barrier). In result, in the region II we obtain: 

+00 

i^lix, t)^ J 9{E - E)e{Vi - E)f3"e^^''-'^*^^dE, ior < x < a. (4) 


Thus the WF in the barrier region constructed by such way, is an analytic continuation of a relevant expression 
for the WF, corresponding to a similar problem with above-barrier energies, where as a stationary expression 
we select the wave exp(iA:2a^), propagated to the right. 

Let's consider the first step further. One can write expressions for the incident and the reflected WP in 
relation to the first boundary as follows 

+00 _ 

i^znc{x,t) = J g{E - E)e{Vi- E)e'-^''-''^^/^'dE, for a; < 0, 

(5) 

iPl^j{x,t) = J g{E-E)e{Vi-E)A%.e-'^''-'^^/^dE, for a; < 0. 


A sum of these expressions represents the complete WF in the region I, which is dependent on a time. Let's 
require, that this WF and its derivative continuously transform into the WF @ and its derivative at point 
X = (we assume, that the weight amplitude g{E — E) differs weakly at transmitting and reflecting of the 
WP in relation to the barrier boundaries). In result, we obtain two equations, in which one can pass from the 
time-dependent WP to the corresponding stationary WF and obtain the unknown coefficients /3'^ and A'j^. 

At the second step we consider the WP, tunneling in the region II and incident upon the second boundary 
of the barrier at point x = a. It transforms into the WP, transmitted through this boundary and propagated 
in the region III, and into the WP, reflected from the boundary and tunneled back in the region II. For a 
determination of these packets one can use Eq. ^ with account ^ , where as the stationary WF we use: 

'fiLcik, x) = v\r{k, x) = /3"e-«^, for < x < a, 

(pl{k,x) = ^^e**'^, forx>a, (6) 

vlefik,x) = Q!°e^^, for < a; < a. 

Here, for forming an expression for the WP reflected from the boundary, we select an increasing part of the 
stationary solution a''exp(^x) only. Imposing a condition of continuity on the time-dependent WF and its 
derivative at point x = a, we obtain 2 new equations, from which wc find the unknowns coefficients Aj, and a". 

At the third step the WP, tunneling in the region II, is incident upon the first boundary of the barrier. Then 
it transforms into the WP, transmitted through this boundary and propagated further in the region I, and into 
the WP, reflected from boundary and tunneled back in the region II. For a determination of these packets one 
can use Eq. ([1]) with account Eq. ([3]) , where as the stationary WF we use: 

^Lcik,x) = ifl^j{k,x), forO<a;<a, 

vlJ{k,x) = A]^e-'^'-', fora;<0, (7) 
ipl^f{k,x) = (3^6'^'', for < a; < a. 

Using a conditions of continuity for the time-dependent WF and its derivative at point x = 0, we obtain the 
unknowns coefficients and (3^. 

Analyzing further possible processes of the transmission (and the reflection) of the WP through the bound- 
aries of the barrier, we come to a deduction, that any of following steps can be reduced to one of 2 considered 
above. For the unknown coefficients a", /3",A^ and A^, used in expressions for the WP, forming in result of 
some internal reflections from the boundaries, one can obtain the recurrence relations: 

Considering the propagation of the WP by such way, we obtain expressions for the WF on each region 
which can be written through series of multiple WP. Using Eq. (JlJ with account Eq. ([3]) , we determine resultant 
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expressions for the incident, transmitted and reflected WP in relation to the barrier, where one can need to use 
following expressions for the stationary WF: 

e^^^j for a; < 0, 

+00 

IJe*''^, for X > a, 



{k,x) 



J2 A^e''''' 



ipref{k,x) 



ri=0 
+ 00 

n=0 



(9) 



— ikx 



for X < 0. 



Now we consider the WP formed in result of sequential n reflections from the boundaries of the barrier and 
incident upon one of these boundaries at point a; = (i = 1) or at point x = a {i = 2). In result, this WP 
transforms into the WP iplj.{x,t), transmitted through boundary with number i, and into the WP ^/^^j{x^t), 
reflected from this boundary. For an independent on x parts of the stationary WF one can write: 



exp(-^a::) 
exp(z/ca:) 

fir 



fief 



exp(iA:x) ' 

vine 

exp{~£,x) ' 

Vine 



exp{—ikx) 
2 

TFT ~ ^ 

exp(p) 



+ 



vine 



exp{ikx) ' 

Vine 



^ exp(-,^a;) ' 



(10) 



vine 



where the sign 



exp{—ikx) exp(^a:) ' exp(— ^a;) exp(^a;) ' 

-" (or "-") corresponds to the WP, tunneling (or propagating) in a positive (or negative) 



a;-direction and incident upon the boundary with number i. Using and R^, one can precisely describe an 
arbitrary WP which has formed in result of n- multiple reflections, if to know a "path" of its propagation along 

the coefficients T-^ and i?.^ can be obtained. 



the barrier. Using the recurrence relations Eq. 

^1 — 



An 

2 



a" ■ 

on+l 



(11) 



Using the recurrence relations, one can find series of coefficients a", /3", and A^. However, these series 
can be calculated easier, using coefficients T^^ and Rf. Analyzing all possible "paths" of the WP propagations 
along the barrier, we receive: 



E 

n=0 

+00 

E^fl 

+00 

E a" 

+00 

E 

n=0 



Fsub 



where 



= Rt 



T,^R^T,-[i+ EiRtRir] = 

71=1 



klD^ 



1+ Y^iRtRiT] = 

+ 00 

1+ E(^2+^r)" 

i=l 



= k 



{e-e)D- 

l±e"2Ca^ 

_ 2mVi 



2fc(i^- fc)e 
2k{i^ + k) 

^ sub 

2ik£,D+, 



(12) 



(13) 



All series E "^"i E E ^^'^ E obtained using the method of multiple internal refiections, coincide 
with the corresponding coefficients a, /?, At and An of the Eq. calculated by a stationary methods [51f?lll5j. 
Using the following substitution 

iC^k2, (14) 

where k2 ~ j- \/2ni{E — Vi) is a wave number for a case of above-barrier energies, expression for the coefficients 
a", /3", Ay and A]f for each step, expressions for the WF for each step, the total Eqs. (fT2| and (fT3| transform into 
the corresponding expressions for a problem of the particle propagation above this barrier. At the transformation 
of the WP and the time-dependent WF one can need to change a sign of argument at ^-function. Besides the 
following property is fulfilled: 



+00 



4-CXD 



= 1. 



(15) 
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2.2 Tunneling of the particle through arbitrary number of rectangular steps with 
arbitreiry heights and widths 

Now let's consider propagation of the particle (from the left) through potential composed of arbitrary number 

of rectangular steps with arbitrary heights and widths. At first, we analyze a case of motion of the particle 
above the barrier. We have the following stationary wave function: 



a2 e""^-" + /?2 e 



aN-i e 
At e^^^, 



ikM-ix _|_ ^-ikN-ix ^ 



at X < xi, 

a.t xi < X < X2, 

a.t X2 < X < X3, 

at xn-2 < a; < Xf^-i, 
at a; > x^-i, 



(region 1) 
(region 2) 
(region 3) 



(region N-1); 
(region N) 



(16) 



where aj and (3j are unknown coefficients, At and An are unknown amplitudes of transmission and reflection. 

The method of multiple internal reflections allows to describe propagation of the particle in such field with 
barriers and to find the transmitted and reflected waves. According to the method, propagation of the particle 
through the barrier is considered by use of wave packet successively by steps of its propagation relatively each 
boundary of the barrier. We conclude that any further step in such consideration will be similar to one from the 
first independent 2N — 1 steps. From the analysis of these steps we combine recurrent relations for calculation 
of amplitudes A^"'^ , 5^"^ , a'"^ and /J^"^ for arbitrary step with number n and determine total amplitudes in 
each region. Further, we shall present results of calculation of amplitudes by approach MIR. 



At first, we determine coefficients , ... T/v-i ^^'^ -^11-^2 • • • -R 



N-i- 



n 



2k 



■3+1 



We find the following recurrent relations: 



'i+l g2ikjXj 



kj+i — k 



R~ — - 



■J ^-2ikj + iXj ^ 



(17) 



R3+1 - ^j+i 



^ m=l 



m=l ' 



R~. 



+ 



1 - Rf+iRj 



- Rj+iRj 



and select the following coefficients as starting: 



R-N-l 



Using them, we calculate successively the coefficients -R^_2 
Now we determine coefficients /Sj: 



R^, f+ = T+. 
. . . , i?2 ■ ■ ■ R 



N- 



and T2 . . . Tj 



/3,- = i;ti (1 + E (^;^7-i)'") = — 

V TO=1 / 1 - 

and amplitudes of transmission and reflection: 



f+ 



Rj R3-1 



Ar = Rt. 



N-1- 



(18) 



(19) 



(20) 



(21) 



After comparison of the amplitudes of wave function obtained by such way with results obtained by standard 
direct method of quantum mechanics for the potential, composed of two rectangular steps of arbitrary height 
we obtain exact coincidence. Increasing of number of steps confirms coincidence of such results. In particular, 
we check a property: 

|At|' + lAnf = 1. (22) 
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All this sure us in that the presented above method MIR gives exact solutions of all amplitudes of wave function 
in all region of its definition. 

Now we consider a case when energy of the particle is less then the height of one step with number m. Then 
for description of propagation of packet inside such step we must use such transformation: 

^ * Cm- (23) 

In this case, we obtain also full coincidence between all amplitudes calculated by method MIR and by standard 
direct method. So, we generalize the method MIR for description of tunneling of the particle through multi-steps 
potential with arbitrary number of steps higher then the particle energy. 



3 Tunneling of the particle through a spherically symmetric barriers 



The one-dimensional formalism of the method of multiple internal reflections in Sec. [5] can be generalized for 
description of a motion of a particle (with its possible tunneling) in spherically symmetric field with barrier, 
allowing to study evolution of nuclear collisions and a-decay of nuclei in spherically symmetric consideration. 



3.1 Scattering of the particle in central field with radial rectangular barrier 



A problem of scattering of the particle on nucleus in spherically symmetric approximation is reduced to the 
problem of scattering of the particle with reduced mass inside the spherically symmetric field with barrier [?] . 
We shall present spherically symmetric formalism of multiple internal reflections for the problem with barrier 
of the simplest form (see Fig. [2]): 



-Vq, at r < Ri; (region I); 

V{r) = { Vi, at i?i < r < i?2, (region II); 

0, at r > i?2, (region III). 



(24) 



We note the region I for < r < region II for i?i < r < R2, and region III for r > R2. We assume that the 

V(r), 



Vi 

d 

-Vo 



I 



II 



III 

Xinc 
-< 

Xtr Xref 



Ri 



Figure 2: Scattering of the particle on the radial rectangular barrier with its tunneling 

particle being under the action of central force V{r) is incident on the external boundary R2 outside. We shall 
consider a case when orbital moment is / = and total energy E is less then the barrier height Vi : E <Vi. 



3.1.1 Transmitted, reflected packets and S-matrix 

At first, we describe standard approach for description of evolution of this process. Solving stationary Schrodinger 
equation in each region, we find wave function: 



), at r < 



(region I), 



at i?i < r < i?2, (region II), 
at r > i?2, (region III) 



(25) 
(26) 
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where Yim.{0,ip) is spherical function, k — ^y/2mE, ki — j^yj2'm{E + Vq), ^ = j-y^2m{Vi — E) (here, we use 
boundary condition of the radial WF at r = to be finite, and its part describing incident of the particle on the 
barrier in region III is normalized on 1). The unknown coefficients S*, A, a and (3 are determined from condition 
of continuity of WF and its derivative at each boundary Ri or i?2 . 

The standard approach to description of evolution of scattering of the particle on barrier consists in use of 
non-stationary wave packets, constructed on the basis of stationary components xQ^tt). Taking into account 
only sub-barrier energies, we obtain: 

-t-oo 

X(r,i)= j g{E-E)e{Vi-E)x{k,r)e-'''"''dE, (27) 



2m 

where the second item in (pS)) represent centrifugal energy, which equals to zero at / = 0, and weight amplitude 
g{E — E) and average energy E are defined like one-dimensional case. 

The problem of a motion of a particle in spherically symmetric field is reduced to the radial problem. The 
radial Schrodinger equation and normalization condition for the radial WF ^-!{r) = x('')/^ is reduced to one- 
dimensional Schrodinger equation with centrifugal potential and the normalization condition for x(^) [I]- We 
describe the particle incident on the external boundary in region III by plane wave exp(— i/cr)? and scattering 
particle by the barrier in region III — by outgoing plane wave S'exp(i/cr), which include both wave, described 
reflection of particle from the barrier, and wave, described transmission of the particle through the barrier when 
it tunnels into region I and then propagates back to region III. The reflected and transmitted waves are formed 
only one component S'exp(jA:r) and it is unclear how to separate them. 

Being non-stationary, the method of multiple internal reflections allows to study the process of scattering of 
this particle inside the spherically symmetric field with the barrier in more details. According to the method, 
the scattering of the particle on the barrier is considered by use of non-stationary WP consequently by steps of 
its propagation relatively each boundary of the barrier (like one-dimensional problem). In result of analysis we 
conclude that any step of further propagations of multiple WPs will be like one from only the first 4 independent 
steps. Analyzing these steps, we combine recurrent relations for calculation of the coefficients A^^\ S'^"\ a*^"^ 
and Z?^"'-' for arbitrary step with number n. 

The total non-stationary WF in each region, taking into account multiple reflections of WPs relatively 
boundaries, can be presented in form of series of these WPs. Analyzing all possible "paths" of propagation of 
packets and using formalism of and i?^, we find sums: 

+ 00 \ 

E S(-^ = -E—T^T+iR^ (1 - i?+i?o ) + Tri?o T+) = 

n— 1 ^ sub 

~ ^1 _ 2ikiRi\ 2i{Ri-R2)-2ikR2 



+ 00 rp-rp- .■,^ ikR.2+lkxR.x-^{R.-2-Rl 



„=o Fsub Fsubik + iO{ki + iS,) 

2^ri + ^e-^-^)e-(«-^-)-^ (29) 

n=0 Fsub Fsubik + 

-t-oc 



n—0 R2 ■^su\ 

2k 



~ ^1 _ g2ifcifli\ 5(2i?,i-fl,2)-ifcfl2 



ii + ki 



Fsub{k + ii) 



where 



Fsub = (1 - i?fi?o )(1 - RtRi) - RtT,-R^T+ = 1 + :^L^e2''=^^i 



(fc i£,){ki "iQ ^-2i(R2-Ri) _ ^ ^g-25(_R2--Ri)+2ifci_Ri 

{k + i^){ki + i^) k + i^ 
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+ 00 

/ 9{E- 




E)Q{Vi 


- E)e-'^'-'^^/^dE, 


Xtr{r,t) = 


+ 00 

/ 9{E- 




E)Q{Vi 


- E)Stre''''-'^^/^dE, 


Xref{r,t) = 


+ 00 

/ 9iE- 





- E)Srefe^^''-^'^^/^dE. 



k + it, + k 

^ i^ + ki ' ^ a(") + fci 

= 0, = 1, (31) 

1 zC + fci ' ^(") z^ + fci 

Here, the coefficients Rf are defined relatively boundary with number z (z = for r = 0, z = 1 for r = i?i 
and i = 2 for r = R2) and are calculated by use of recurrent expressions between 5^"), a^") and 

Now we define incident, transmitted and reflected WPs relatively the barrier (in region III): 



(32) 



where 

+00 

S=Str + Sref. Str = S^-^\ Sref = S^°\ (33) 

71=1 

So, method MIR separates S- matrix into two components, corresponding to the amplitudes of the transmitted 
and reflected WPs relatively the barrier. This property fulfills for any energy level and allows to find coefficients 
of penetrability and reflection of the particle relatively the barrier (and also probabilities of processes of scattering 
through compound nucleus formation and potential scattering). 

The coefBcients S^'^\ y4("\ a^") and for each step, WFs for each step, the coefficients and Rf, total 
sums of the coefficients S'("\ A^"'\ a^") and /J^"-* at ^4]) transform into corresponding expressions for motion 
of the particle above the barrier. The total sums ([29]) of the coefficients S'^"'\ ^("\ a^") and coincide with 
corresponding coefficients S, A, a and (3, calculated in the direct method of QM. We have: 

\S\ = 1. (34) 

3.1.2 Phase times of tunneling and reflection relatively the barrier 

According to method of stationary phase, we have 

d d d 

^arg Xinc{r, t) ^ ^arg xtr(r, t) = ^arg XTef{r, t) = const. (35) 



Now we consider the first step. Let maximum of WP in region III be incident to the external boundary of the 

^ref 



barrier at point i?2 at time moment tine- Using (|35|) . we define time moment t^^X of leaving of the maximum of 



the reflected WP outside in region III: 

(1) , 2mR2 , ^ aarg^W 

For time moment il"-* of leaving from the barrier region outside of the maximum of rt-multiple WP we write: 

(„) , 2mR2 , ^ aarg5(") 

+ ^ + (37) 

Using ([35)1 at point r ~ R2 and taking into account ((36| , (IST)) , we find time of transmission of WP through the 
barrier and time of reflection of WP from the barrier: 

_ _ _ 2mR2 (9arg Str 

Ttr-ttr t,nc- + ^ : 

_ _ _ 2mR2 , J&TgSref ^'^^^ 

Tref — tref tine — + ^ qj^ 
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For WP tunneling under the barrier we obtain: 



^tun 
Tref 



2m 



(39) 



For propagated WP above the barrier we obtain: 

_ 2m(i?2 - i?i) 



' prop 



^d_ ^ fc2 - fci - (fc2 + fci)e^'^i^^ 

We"^^^ {k + k2){ki + k2)Fabove ' 



(40) 



Tref — 0, 

where Fatove can be obtained from Fsub at (1141) . Also we have 



Ttr = 
Tref 



hk 

2mR2 
hk 



darg S d&Tg 7 
_ ^ d'Ai-g S ^ ^9arg (1 - 7) 



dE 



dE 



where 



7 



Str 



(41) 



(42) 



One can see that tunneling and reflection times are different only by one item [?] . 

Now we consider scattering of the particle on enough high and wide barrier. In such case for sub-barrier 
tunneling we obtain: 

_ 2m 4mi?isin2fcii?i(l - 2cos2A;ii?i) 

" ^ + ae(i-cos2fcii?i) • ^ ' 



Here, the following consequence of approximations are used: ^(i?2 — Ri) 
3.2 a-decay 



+00, ^ +00, i?2 — Rl 



Now let's generalize the method MIR for description of a-decay in spherically symmetric approximation. We 
shall assume that starting from some time moment (one can name it as time moment ^formation of formation 
of the a-particle and daughter nucleus) the decaying nucleus can be considered as composite system of the a- 
particle and daughter nucleus. We come from a problem of motion of the a-particle in the field of the daughter 
nucleus to another problem of propagation of the particle with reduced mass m in spherical symmetric field 
with a radial barrier. 

For simplicity, we shall consider a case when total energy of system E is higher then the barrier height Vi : 
E > Vi. For forming a main idea of application of the multiple internal reflections for description of a-decay, 
we shall assume that the radial potential V(r) can be separated into 3 regions: internal region I, region of the 
barrier II, and external region III, and in each region we know exact general solution of the stationary wave 

V(r), 



Vi 



t ^ 


II 


III 

Xtr 

► 


Zref 




-« 

Xinc 




R] R2 r 







Figure 3: Propagation of wave packet from the internal region outside with its tunneling through the barrier in the 
a-decay 



function in the form: 



(44) 
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Aic^ (r) + Bic+(r), at < r < i?i (region I), 
X{r) = { A2C-{r) + B2C+{r), at i?i < r < i?2 (region II), (45) 
B3C^{r), at r > i?2 (region III) 

where c~(r) and cf{r) are ingoing and outgoing waves, Ai and Bi are constant coefficients. According to 
ideology of multiple internal reflections, we describe the a-decay by use of wave packet describing the particle 
with reduced mass, successively by steps of its propagation from the internal region I outside through the barrier 
region II relatively each boundary (see Fig.[3|). 

In the first step, we consider beginning of the a-decay as propagation of the packet in the region I outside 
at t > tformation, which at somc time moment tine is incident on the internal boundary of the barrier at We 
name this time moment as time moment of beginning of the a-decay. The incident WP forms new WP in the 
region II at time moment ttr propagating from the boundary i?i outside, and new WP in the region I at time 
moment tref propagating from this boundary inside back. We define the incident, transmitted and reflected 
packets concerning the boundary Ri so: 

X^nc{r, t) ^ J g{E - E)Q{E - Vi)x;nc{r, k)e-'^*^^dE, 



xi'Jir^t) = J g{E~E)eiE-V,)x[lHr,k)e-'^'/''dE, (46) 



-\-oo _ 

xl'j/(r, t) = J g{E - E)Q{E - V^i)xlt/(^> k)e-'''"''dE 



where g{E — E) is gaussian for above-barrier energies and we have: 

X^„c (r, k) = b[^^c+ (r) , at < r < i?i , 

xil^ {r,k)^ Bi'^c+ (r) , at i?i < r < i?2 , (47) 
Xrefir, k) = A^pc^ir), at < r < i?i 

Also we shall neglect by variation of weight amplitude 

g_(£;)^g«(i?)^5S(^)- (48) 

Using continuity condition for total wave function composed on these WPs and its derivative at Ri , we find 
the unknown coefficients A^^'^ and B^2^: 

_ (,) C+(r)|-4(r)-4(r)^c+(r) 

1 ~ 1 d _ _ d ' 

4W^cr(r)-cr(r)-c+(.) ^^^^ 

.(1) _ pd) ^^^"^^;'^^^'^)'^^^'^)^^^^^"^ 

2 ~ ^1 Q _ _ d 

4W^Ci (r)-Ci (r)— c+(r) 

Then, like one-dimensional problem (see Sec. 12. ip and the scattering problem of the particle on nucleus (see 
Sec. 13. ip . we consider further propagation of WP through the barrier outside in next steps. Here, any nex 
step with propagation and reflection of WPs concerning the barrier boundaries can be reduced to one from 
4 independent steps. For determination of the unknown and one can find recurrent relations and 
calculate the coefficients and Rf: 
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B. 



(1) 



(1) 



(1) 



Rn 



0, 

-1, 



T+ 



ctiR2)g^C^{R2) - C2-(i?2)^C+(i?2) 

C+(i?2)|:C2-(i?2) - C2-(i?2)|:C+(i?2)' 
4(^2)^C3+(j?2)- 4(^2)^4(^2) 

4(^2)|-C2-(i?2) - C2-(i?2)|-C+(i?2) 



c+iRi 



dr 



c-{R,)-c-{R,)—4iR,) 



(50) 



4{Ri)^c^{Ri)-c^{Ri)-^c+{R^) 
c^{Ri)—c^iR,) ~ c^iRi)—c^iRi) 



d 



d 



4{Ri)—c^{R,) - c^{R,)—4{R,) 



dr dr 
One can find total transmitted and reflected WPs relatively the barrier: 



+00 





+ OC 



Xtrir, t) 

Xrefir, t) 

All sums of the coefficients A]"' and B 



J g{E - E)Q{E - Vi) YZ=i b'^^'^ c+{r)e'^'-'^'/^dE 





(") n(" 



/ g{E-E)Q{E-V,)Y.l^,A'C'<^-A^ 



*/^dE. 



(51) 



n=l 



n=2 

+ 00 , 

n=l 



(n 



E ^2 

+00 



(n 



are calculated on the basis of and . We obtain: 



-B\ 
= B\ 



Rq (1 - RiRq) 

1 



1 — Ri Rq 



R^{1^R+Rq){1-RYR+) 



(52) 



n=l 



where 



(i-i?+i?o-)(i-i?ri?j)' 

il-R+RQ){l-R^R+) 



D+ _ p+ I ^1 ''l ^2 n+ _ n- 

1 — Ri R2 1 — Ri Rq 



(53) 



All expressions for sums ^ A^"'' and ^ B^"'' obtained by method MIR coincide with the corresponding coef- 
ficients Ai and Bi calculated by standard direct method (see also [?]). If the total energy E is less then the 
barrier height Vi , then expressions for sums of the coefficients A^"^ and B^"'' are obtained by use of ([49l) , (l50l) , 



(1521) . ([53|) and (fH)) for the barrier region II. These coefficients for arbitrary step are calculated on the basis of 
(j49| and the coefficients and Rf. Amplitude b[^^ is found from the normalization condition for the radial 
WF. 

A non-stationary analysis of the a-decay can be fulfilled, like one-dimensional problem (see [SI llOj ) or 
the scattering problem (see Sec. 13. ip . We shall estimate duration of the a-decay on the basis of the barrier 
penetrability. In consideration of motion of the a-particle from the internal region outside with its tunneling it 
needs to take into account multiple internal reflections inside the internal region I, formed in result of reflections 
of packets from the barrier inside. But, the barrier penetrability can be related to duration of decay and to 
multiple reflections inside the internal region by different ways. Therefore, we give two independent approaches 
for deflnition of the penetrability coefflcient T of the barrier, describing propagation of the particle from the 
internal region I into external region III. 



Approach A: 
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In the first case, we define the coefficients of penetrability and reflection with taking into account of multi- 
ple internal reflections of packets inside the internal region. Here, as the incident, transmitted and reflected 
stationary parts of we use Xincir) from (|T7l) and Xtr{r), Xref{r) from ([?T|) . We have: 



T 



R 



+ 00 

11=2 



(n) 
3 



(1) 



+ 00 

n=l 



in) 



c+(r)|:c+(r)* 



4irr§;ct{r) 



(54) 



(1) 



d 



d 



ct{r)^c+{rr ~ct{rr-c+{r) 



dr 



At such definition, the coefficients T and R are not separated evidently into factor, characterized the penetra- 
bility of the barrier directly without taking into account multiple reflections of packets inside the internal region 
I, and factor characterized such multiple internal reflections. 

Here, we introduce duration of a- decay Tdecay, defining it as difference between time moment ttr of leaving 
of the total transmitted WP from the barrier (time moment of finishing of a-decay) and time moment tine 
of incident of WP inside region I on the internal boundary i?i in the first step (time moment of beginning of 
a-decay). Also, we define time moment of reflection of the particle from the barrier Tref as difference between 
time moment t^ef of reflection of the total reflected WP from the barrier into the internal region I and tine- 
Using (|35|). we obtain: 



Ttr 



Tref 



tt. 



dE 

tref 

d_ 
'dE 



h 

+ 00 



_d_ 

dE 



arg c+(i?2) 



_d_ 

dE 



arg c^(i?i)- 



tine — h 



dE 



d _ d 
—arg cr(i?i) - ^arg c+(i?i) 



(55) 



arg 



n—l 



(n) 



^arg B[ 



Approach B: 

In the second case, we define the penetrability coefficient of the barrier without taking into account of multiple 
internal reflections of packets inside the internal region, i. e. as in one-dimensional case. However, such internal 
reflections we include inside the internal region, like in the semiclassical way of definition of half-lives — through 
introduction of independent normalized factor F defined a number of "collisions" of the a-particle with the 
barrier inside the internal region I. As non-stationary characteristic of the a-decay, we define half-life r on the 
basis of width T of the a-decay. 

3.2.1 Decay with radial a-nucleus potential in a form of one rectangular barrier and rectangular 
well in the internal region 

As an example of application of the method MIR to calculation of the barrier penetrability and estimation of 
time characteristics in approach A, we shall consider the a-decay with the barrier of the simplest form: 

{~Vo, at r < Ri; (region I); 

Vi, ati?i<r<i?2, (region II); (56) 
0, at r > i?2, (region III). 

Let's orbital moment is Z = and total energy E less then the barrier height. Stationary WF is: 

i;{r,9,^)^^YU0,^), (57) 
r 

^(g-ifeir _ gifcir)^ atr < Ri, (region I), 

Xir) = { ae«'' + /3e~€^ at Ri < r < R2, (region II), (58) 
S'e^'=^ at r > i?2, (region III) 
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where Yi„i{9,ip) is spherical function, k — ^\/2mE, ki = ^y/2m{E + Vq), ^ — j-y^2m{Vi — E). Taking into 
account we write: 

cf = e±*'=l^ c+ = e'''''. (59) 
The coefficients T^^ and Rf are calculated on the basis of (|50p and using form of cf: 

, _ 2kl fR^+ikiRi p+ _ ^1 ~ ^^ 2lkiRi 

It, + ki ki+ It, 

. _ 2t^ ^^R^_ikB.2 R+ _ !i_^„-2efl2 

2_"fc + iC ' l~i^ + k ' (60) 

Tg =0, = 

7^- ^ ^eRi+ikiRi ^- ^ »g - fcl ^2;fli 

^ iC + ^i ' ^ iC + fci 

The coefficients of penetrability T and reflection i? of the particle relatively the barrier can be obtained from 
([54l) . Time of tunneling ( ) and time of reflection can be found from ([55]) . We obtain: 

d d 
Ttr = ^-Q^ arg Ctr, Tref = arg C^e/ ■ (61) 



where 



iAki^ 



(fc + iC)(^i + + (fcl - iS,)exp{2ikiRi)) ' 



_ (^^ - k){ii - fcl - (z^ + ki)exp{2ikiRi)) _ 
ii^ + kM + k,-izi-k^)expi2zk,R,)f''^^'^^^' 



ref 



_1_ 



- |^e.p(2.fci^0 ^ ^ ~ ^ - g e..(2^(i^r - R.))) + (62) 
4zfce(fc-*0 -exp(2e(i?i-i?2) + 2zfcii?i); 



{t^ + k)^i^ + ki) 




C2 = 1 + '^^exp{2^kR,) 1 - " f !S^^, " f ea;p(2^(i?i - i?^)) 



fci + i^ («c + fc)(ie+fci) 

4ifc^(fc - i^) 



(i^ + k)^^^ + ki) 



exp{2^{Ri - R2) + 2ikiRi). 



3.2.2 Decay with the radial a-nucleus potential composed of arbitrary number of rectangular 
steps 

Now we shall consider the a-decay with a realistic a-nucleus radial barrier, when this barrier is enough well 
approximated by number of rectangular potential steps. The first region I we define starting from point i^min, 
and we assume that the a-particle is formed here and moves outside. Supposing that a-particle can be formed 
enough far from zero, we use non-zero value for i?min- As before, we shall analyze the motion of this particle on 
the basis of multiple internal reflections of WP successively by steps of its propagation relatively all boundaries 
between steps in the multi-step radial barrier. In consideration of multiple internal reflections of multiple packets 
inside radial internal region starting from zero and in dynamical estimation of the a-decay, we use the second 
approach B and shall find half-live of the a-decay. Then comparing obtained results with half-lives obtained by 
the semiclassical approach, we shall estimate how the method MIR can be effective. 
At first, we write total wave function: 

e*'"' -I- Ar e"*'"'', at i?,ni„ <r <ri, (region 1) 

a2 e*'^^'' -I- e"*'^^'', at ri < r < r2, (region 2) 

as e*''-'"" + (3^ e"**^^'', at 7-2 < r < 7-3, (region 3) 



(63) 



a„_le*'="-l'■ + /3^^_le-*'="-l^ a.t rN-2 < r < r^-i, (region N-1); 
At at VN-i <r < i?„iax, (region N) 



where aj and Pj are unknown amplitudes, At and Ar are unknown amplitudes of transmission and reflection. 
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Like one-dimensional problem, we define coefficients , . . . T^_i and , i?J . . . R'^_i- 

% + fcj + i ' /Cj + fcj + i 

fcj + fcj + i ' fcj + fcj + i 

Using recurrent relations: 



(64) 



^7+1 = ^7+1 + t;>i^7^>i (i + E (^;+i^7)"0 = ^7+1 + 7%rg"- ' (65) 

^ m=l ' i — nj_^_iKj 

^ m=l ^ t — ^j+l^j 

and selecting as starting the following values: 

^jv-i = -fi'jv-i) -^1 ^ Ri ^ = ^ (66) 

we calculate successively coefficients Rjf_2 ■ ■ ■ Rt ^ ^7 • • ■ -^jv-i ^^"^ -^2^ • ■ ■ -^jv-i- ffiiishing, we determine 
coefficients Pj: 

the amplitudes of transmission and reflection: 

AT^f+_^, Aji^R+ (68) 
and corresponding coefficients of penetrability T and reflection R: 

Tmir = 7— \ At\ , Rmir ~ l^fll • (69) 
As test of the found solutions, we use the property: 

^ |^T|' + |^i?|' = l or Tmir + Rmir = 1- (70) 
ki 

3.2.3 Width T and half-live 

We define width T of the a-decay, in semiclassical approximation, by following the procedure of Gurvitz and 
Kdlhermann [21j : 

r = P„ F — T 71 

4m 

where Pa is the a-particle preformation probability and F is the normalization factor. T is the penetrability 
coefficient in propagation of the particle from the internal region outside with its tunneling through the barrier, 
which we shall calculate by approach MIR ar by approach WKB. In approach MIR this coefficient we define so: 



TwKB = cxp < -2 
R2 



(Q-V{r)) dr\ (72) 



where R2 and i?3 are the second and third turning points. According to [13', the normalization factor F is 
given by simplified way so: 

f R-2 \ -1 
dr 
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or by improved way so: 



I J li(r) 

^Ri 

The a-decay half-live is related to the width T by well known expression 




k{r') dr' - J 



dr \ . (74) 



Ti/2 - ?i ln2/r. (75) 



4 a-nucleus potential 



For description of interaction between a-particle and the daughter nucleus we shall use the a-nucleus potential 
in the general form: 

V{r, 0, 1, Q) = vc{r, 0) + VN{r, 0, Q) + vi{r) (76) 

where vc{r,0), VN{r,0,Q) and vi{r) are Coulomb, nuclear and centrifugal components. There are different 
approaches in definition of parameters for these components. In this paper we shall use two variants. 

4.1 a-nucleus potential of V. Yu. Denisov and H. Ikezoe 

In the first case, we shall use the approach proposed by V. Yu. Denisov and H. Ikezoe in [IB] (for simplicity, we 
shall name the potential (|76p with parameters defined according to [1.6j as potential in the form DI). According 
to (6)-(10) in [16j, Coulomb vc{r,0), nuclear VN{r,0,Q) and centrifugal vi{r) components have the following 
form: 

( 1 + -^P2Y2o{0) J , for r > r™. 



2Ze^ ( 3 H 3R- 
VNir,0,Q) = — 



3 N 

■P2Y2o{e)(2~^)\, forr<r 
ViA,Z,Q) 



r-r„(0)' (78) 

1 + exp 

a 

We define the parameters of the Coulomb and nuclear components as (see relations (14), (16)-(19) in |16j): 

V{A,Z,Q) = -(30.275 - 0.45838^/^1/=' + 58.270/- 0.24244 Q), (80) 

R = i?p(l + 3.0909/i?2) + 0.1243 (81) 

Rp = 1.24^1/3 (1 + 1.646/A- 0.191/), (82) 

t = /-0.4yl/(yl + 200), (83) 

d = 0.49290, (84) 

/ {A-2Z)/A. (85) 
According to relations (21)-(22) in [T^, we also use: 

r„,(0) = 1.5268 + /?(0), Ri0) = R (1 + f32Y2oi0)). (86) 

Here, A and Z are the nucleon and proton numbers of the daughter nucleus, respectively; Q is the Q-value, for 
the a-decay, R is the radius of the daughter nucleus, V{A, Z, Q, 0) is the strength of the nuclear component; 
Tm is the effective radius of the nuclear component, d is the parameter of the diffuseness; l2o(^) is the spherical 
harmonic function of the second order, is the angle between the direction of the leaving a-particle and the 
axis of the axial symmetry of the daughter nucleus; (32 is the parameter of the quadruple deformation of the 
daughter nucleus. 
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In this paper we restrict ourselves by spherical symmetric approximation of a-decay. Then the Coulomb 
and nucleus components are transformed into the following: 

, for r > r„, 

vc{r,9) = { r ^2. (87) 
S 3 5- >, for r < Tm, 

^1^4^ ,88, 



1 + exp ■ 



d 

where 

r„ = 1.5268 + (89) 

4.2 a-nucleus potential of B. Buck, A. C. Merchant and S. P. Perez 

In the second case, we shall use approach proposed by B. Buck, A. C. Merchant and S. P. Perez in JTT (we 
shall name the potential (|76p with the parameters defined according to ^13j as potential in the form BMP). 
According to (l)-(3) in [13], the nuclear component is: 

, . I + cosh Rbmp /d 

VN{r) = -Vo , , , , 7-^ Ti, 90 

coshr/d + cosh RBAip/d 

Coulomb component is: 

2Ze2 



Voir) = < 

and centrifugal component is: 



for r > Rbmp, 

ZeV r2 \ (91) 



3 - I , for r < Rbmp 

BMP . 



{I + i)2 



- (92) 

Values Vq, Rbmp and d for the different nuclei are presented in tables in [13] . 



5 Analysis of the method MIR in a-decay 

We shall study, how the presented above method of multiple internal reflections really works (how much it is 
effective, accurate, gives stable result) in the problem of determination of half-lives in a-decay of heavy nuclei, 
which can be considered in enough good approximation as spherically symmetric. Today, there is a lot of well 
developed methods of calculations of half-lives of a-decay, they are experimentally studied well. So, for method 
MIR analysis we have rich theoretical and experimental material in this problem. 

We shall use two a-active nuclei: ^^"Po and ^^*Po. Such a choice we explain by they have small coefficient 
of quadruple deformation (32 and in enough good approximation can be considered as spherical. They have 
very similar shapes of barriers and parameters of a-decay, but their half-lives (and Q-values) are different 
essentially (according to jl3 : roxp(^^"Po) = 1,2 • lO'^ sec and Tcxp(^^'*Po) = 1,6 • 10"^ sec). We shah be 
interesting in how accurately the method MIR allows to calculate these values and to describe such difference 
(in a comparison with the semiclassical approach). Besides, we presently have numerical methods of calculations 
of wave functions with high accuracy tested for these nuclei (that is shown essentially in calculations of spectra 
of photon bremsstrahlung accompanying a-decay, see fl71 118 1 ^19 ] and computer programs. 

At first, we shall consider a-decay of ^^^Po. We shall use a-nucleus potential in the spherically symmetric 
BMP form. We shall study a-decay on the basis of leaving of the particle with reduced mass from the internal 
region outside with its tunneling through the barrier inside the region of r with internal boundary -Rmin (we 
assume it to be very close to zero) and the external boundary i?niax (we assume it to be very large). To apply 
the method MIR for description of this process, we change the radial a-nucleus potential inside the studied 
region into approximated one, which consists from finite number N of rectangular potential steps. We define 
the region of formation of the particle inside the first interval with i?niin ^ t < Ri (here we assume: Q > V(r))^ 
from here, this particle starts to move outside. According to method MIR, we consider leaving of the particle 
on the basis of WP consequently by steps of its propagation relatively each boundary. Using technique of the 
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Number N 


Transmission amplitude 


Reflection amplitude 


Penetrability 


Half-live 


of intervals 




5 (At) 




^{Ar) 


T 


At 


ri/2, sec 


1000 


6,97E-10 


-l,200E-9 


-0,976 


0,213 


1,2946E-18 


- 


l,55E-4 


2000 


6,612E-10 


-l,148E-9 


-0,9683 


0,249 


1,1790E-18 


0,1156E-18 


l,62E-4 


3000 


6,504E-10 


-l,140E-9 


-0,9636 


0,267 


1,1580E-18 


0,0210E-18 


l,87E-4 


4000 


6,452E-10 


-l,139E-9 


-0,9608 


0,277 


1,1507E-18 


0,0073E-18 


l,94E-4 


5000 


6,422E-10 


-l,1384E-9 


-0,9590 


0,2832 


1,1473E-18 


0,0034E-18 


l,77E-4 


6000 


6,403E-10 


-l,13834E-9 


-0,9577 


0,2874 


1,1454E-18 


0,0019E-18 


l,78E-4 


7000 


6,389E-10 


-l,13838E-9 


-0,9568 


0,2905 


1,1443E-18 


0,0011E-18 


l,72E-4 


8000 


6,379E-10 


-l,13848E-9 


-0,95613 


0,2929 


1,1436E-18 


0,0007E-18 


l,66E-4 


9000 


6,371E-10 


-l,13859E-9 


-0,95556 


0,2947 


1,1431E-18 


0,0005E-18 


l,77E-4 


10000 


6,365E-10 


-l,13871E-9 


-0,95510 


0,2962 


1,1428E-18 


0,0003E-18 


l,79E-4 



Table 1: Dependence of the amplitudes of transmission and reflection, penetrability coefficient and half- live on 
number N of intervals in the region from i?min up to i?max for a-decay of ^^''Po (we use: i?min — 0, 11 fm, 
-Rmax = 100 fm; At = |7n+i — Tn\ is parameter of convergence in calculations of the penetrability coefficient T 
where T„ and T„+i are the penetrability coefficients for each previous and further numbers N; accuracy of the 
obtained coefficients is |1 — Tn — Rn\ < 1, 5 • 10^^^ for each value N, that gives the first 14-15 digits as rehable 
for amplitudes and coefhcients) 

coefficients and Rf in (jMll -iP l) . we calculate amplitudes Pj in each interval by ([S7|) . total amplitudes of 
transmission At and reflection Ar by (|68p . The penetrability coefficient Tmir, describing leaving of the particle 
from the internal region outside with its tunneling in the region from i?niin up to i?max) we calculate by ([55]) . 
In this paper we restrict ourselves by definition ([75]) for the normalization coefficient F2 in WKB-approach, we 
find width F by ([TTjl and half- live tmir by ((75)l . 

Further, we define the penetrability coefficient Twkb in the WKB-approach by (|72p . on the basis of which 
we recalculate F- width and half- live twkb by (|71l) and (TTS]) . Comparing the found half- lives (also amplitudes 
of transmission and reflection, penetrability coefficients) by two approaches, we analyze the effectiveness of the 
method MIR for the studied nucleus. 

5.1 Analysis of convergence of the method MIR in dependence on the number of 
intervals inside the studied region 

Let's analyze whether the method MIR gives unique values for the penetrability coefficient, amplitudes of 
transmission (in the last interval with number N) and reflection (in the first interval with number 1) after taking 
into account all needed multiple internal reflections relatively all boundaries, in dependence on the number N of 
intervals inside the region from i?min to i?max- According to [13], for a-decay of ^^^Po we use Qa = 7, 865 MeV, 
-Rbmp = 7, 417 fm. The calculated penetrability coefficient T with increasing of the number N is presented in 
table [H Here, one can see that the penetrability coefficient T tends to 1, 142 • 10~^^ which can be considered as 
limit value. Parameter At = |Tn-i-i — Tn\ demonstrates a convergence of such calculations enough well, which 
is decreased with increasing of N. We include into the table also the amplitudes of transmission and reflection, 
which have the flrst some digits stable. We establish (at first time): 

• The method MIR has a convergent algorithm of determination of the amplitudes of wave function, the 
coefficients of penetrability and reflection which are unique and can be considered as reliable (inside 
interesting accuracy) concerning to the realistic a-nucleus potential in the form BMP. 

• The method MIR determines the half-live for •^^"'Po in the interval tmir = 1, 55 • 10^'' — 1, 79 • 10"^ sec 
which is essentially closer to experimental value Toxp = 1,6- 10^^ sec in a comparison with half-life in the 
approach WKB, which equals twkb = 1,1- 10"** sec. 

5.2 How is the barrier penetrability changed at increasing of the external bound- 
ary -Rmax? 

Let's analyze whether the convergence in determination of the penetrability coefficient by approach MIR at 
increasing of the external boundary -Rmax exists. In approach WKB the region of the barrier located between two 
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Boundary 


Transmission amplitude 


Reflection amplitude 


Penetrability 


Half-live 






5 (At) 




^(Ar) 


T 


At 


Ti/2, sec 


31 


l,15E-9 


-3,27E-9 


-0,95509 


0,29630 


1,6970E-18 




l,17E-4 


40 


-7,94E-10 


-l,47E-9 


-0,95509 


0,29629 


1,1258E-18 


0,5712E-18 


l,88E-4 


50 


-1,47E-10 


-2,63E-9 


-0,95509 


0,29629 


1,1326E-18 


0,0068E-18 


l,75E-4 


60 


-9,51E-10 


-l,41E-9 


-0,95509 


0,29628 


1,1460E-18 


0,0066E-18 


l,73E-4 


70 


9,56E-10 


9,87E-10 


-0,95509 


0,296282 


1,1463E-18 


0,0003E-18 


l,85E-4 


80 


-9,41E-10 


-9,58E-9 


-0,95510 


0,296282 


1,1452E-18 


0,0011E-18 


l,83E-4 


90 


4,00E-10 


l,25E-9 


-0,95510 


0,296282 


1,1443E-18 


0,0009E-18 


l,70E-4 


100 


6,36E-10 


-l,13E-9 


-0,95510 


0,296282 


1,1428E-18 


0,0015E-18 


l,79E-4 



Table 2: Dependence of the amplitudes of transmission and reflection, penetrability coefficient and half- live on 
the external boundary R^ax for a-decay of ^^''Po (we use: i?min = 0, 11 fm, width of each interval is 0,01 fni 
(100 intervals inside 1 fm); At = |7ri-i-i — Tn\ is parameter of convergence in calculations of the penetrability 
coefficient T where T„ and Tn+i are the penetrability coefficients for each previous and further values of i?max; 
accuracy of the obtained coefficients is |1 — r„ — < 1, 5 • 10"^'"' for each value of i?max) 

turning points i?2 and i?3 makes up a main part in determination of the penetrability (we obtain: i?2 = 8, 18 fm 
and i?3 = 10, 09 fm). Assuming that approach WKB gives enough well results, we shall suppose that in approach 
MIR the main region forming the penetrability is the barrier region, while the internal and external parts of 
the potential do not influence on it practically. Keeping width of each interval to be constant, we shall increase 
gradually the external boundary -Rmax (by increasing number N of intervals) starting from value near to the 
external turning point i?3 (-Rmax > R3) and we shall calculate the amplitudes and the penetrability coefficient. 
These calculated values with different i^max are presented in table [21 One can see that method MIR gives the 
convergent value for the penetrability coefficient at increasing of -Rmax- But such convergence is not so stable 
as in increasing of the number of intervals at fixed -Rmax- However, starting from -Rmax = 60 fm, we obtain the 
first 3 digits stable. Convergence of the amplitude of refiection is higher essentially. We conclude: 

• The amplitude of transmission and the the penetrability coefficient are convergent with increasing of -Rmax 
but not so strong, that can be explained by that not only the barrier region but also the external part of 
the potential up to 60 fm takes place in determination of these characteristics (in contrast with approach 
WKB). 

• Half-live obtained by method MIR inside only the barrier region is ttviir = 1, 17 • lO^'' sec, which is close 
to value TwKB = 1,1' 10^'* sec obtained by approach WKB in [13]. This confirms the conclusion about 
essential influence of the external part of the potential to the penetrability coefficient determination. 

5.3 Dependence of the penetrability on the starting point 

Now we shall analyze how a location of the internal boundary -Rmin has influence on the penetrability. In the 
Fig. [Hone can see that half-live of a-decay of ^^^Po is changed in result of displacement of -Rmin- Taking into 
account that width of each interval is 0,01 fm, we shall consider point -Rmin as a starting point (with error up 
to 0,01 fm), from here the a-particle begins to move outside (to the right) and is incident to internal side of the 
barrier starting the first stage of the a-decay. Analyzing figures, we establish a (non- constant) dependence of 
the penetrability coefficient and half-live on the location of the starting point -Rstait, where the particle starts to 
move outside by approach MIR (it has obtained at the first time). 

This dependence opens a new way of determination of the most probable radial coordinate of formation of 
the a-particle before its motion outside in the first stage of the a-decay on the basis of a comparable analysis 
of the calculated half-live with its experimental value: it needs to find such starting point, when the calculated 
half-live by approach MIR is maximally close to its experimental value. One can assume, that from such space 
point the a-particlc begins to move outside in the first stage of the a-decay, and there the a-particle was formed. 
We give such definition: 

The radial coordinate -Rform of the most probable formation of the a-particle before a-decay is such starting 
point, at which the calculated half-live tmir is closed maximally to its experimental value Toxp. 

For a-decay of ■^^'^Po we obtain -Rform = 8, 02 fm (at -R2 = 8, 18 fm). At such starting point the calculated value 
of half-live by approach MIR is tmir = 1,5786 • 10^'' sec, that is essentially closer to the experimental value 
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Figure 4: Dependence of half-live for "^^''Po on location of the starting point -Rstart: (a) — dependence of half-live on 
the starting point from up to 9 fm; (b) — increased region where half-live is changed slowly; (c) — determination of 
coordinate of formation of a-particle before a-decay, as such starting point Tistart, where the calculated half-live is close 
maximally to its experimental value: we obtain -Rstart ~ 8, 02 fm at tmir ~ 1, 5786 ■ 10"* sec 



Texp = 1)6- 10^'* sec for this nucleus, in a comparison "with half-live twkb = 1,1' 10^^ sec obtained in [13j by 
approach WKB. 

So, on the basis of analysis of a-decay of ^^*Po we have seen that the coordinate of formation of the a- 
particle is located inside internal region and is very close to the second turning point R2- But "whether is this 
result accidental, obtained in result of accidental coincidence of parameters of the potential and energy Ql Let's 
analyze another nucleus — ^^°Po. While these t"wo nuclei have practically similar shapes of their potentials, their 
(experimental) half-lives are different essentially. We shall clarify "whether our algorithm, used in the method 
of multiple internal reflections, allo"ws to determine the coordinate of formation of the a-particle for ^^°Po and 
"where it "will be. In Fig.[5]one can see that method MIR determines this coordinate simply also and it locates like 
nucleus ^^^Po — close to the second turning point R2. We obtain value -Rform = 7, 53 fm (at R2 = 7, 642 fm). 
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Figure 5: Dependence of half-live for ^^"Po on location of the starting point .Rstart: (a) — dependence of half-live on 
the starting point from up to 8 fm; (b) — increased region where half-live is changed slowly; (c) — determination of 
the coordinate of formation of the a-particle before a-decay, as such starting point J?start , where the calculated half-live 
is close maximally to its experimental value: we obtain 7?start ~ 7, 53 fm at tmir ~ 1, 16313 • 10^^ sec 

At such starting point the calculated value of half-live by approach MIR is tmir = 1,1631 • 10+^ sec, that 
that is essentially closer to the experimental value Texp = 1,2 • 10"'"'' sec for this nucleus, in a comparison "with 
half-live twkb = 3,6- 10"*"^ sec obtained in [T3] by approach WKB. So, the found conclusion about formation 
of the a-particle (and also improvement of agreement bet"ween the calculated and experimental half-lives after 
consideration of leaving of the a-particle from the coordinate of its formation by approach MIR) is stable both 
for short-lived and for long-lived nuclei. 

5.4 Results of calculation of the coordinate of formation of the a-particle for nuclei 
with charge number Z = 84 

The calculated coordinates i?form of formation of the a-particle for the a-active nuclei "with charge number 84 
are presented in table [3l One can see that the coordinate -Rform is located at small distance Ar, near the second 
turning point i?2. At such choice of the starting point the calculated half-live tmir by approach MIR for each 
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Parent nucleus 


Half- live- values , 


sec 


Turnin 


g points 


a-particle formation 


A 


Q, MeV 


-Rbmp, fm 


TWKB 


T MIR 




i?2, fm 


i?3, fm 


-Rform, fm 


Ar, fm 


201 


5,830 


6,956 


6,6 E4 


5,7101 E4 


5,7E4 


7,6350 


40,4543 


7,5905 


0,0445 


202 


5,733 


6,959 


1,2 E5 


1.3016 E5 


1,3E5 


7,6363 


41,1803 


7,5760 


0,0603 


203 


5,524 


6,965 


2,1 E6 


1.9000 E6 


1,9E6 


7,6385 


42,7536 


7,5972 


0,0413 


204 


5,516 


6,965 


1,4E6 


1,9028 E6 


1,9E6 


7,6384 


42,8182 


7,5911 


0,0473 


205 


5,356 


6,970 


1,7E7 


1,0003 E7 


1,0E7 


7,6405 


44,0971 


7,5723 


0,0682 


206 


5,358 


6,970 


9,8 E6 


1,4036 E7 


1,4E7 


7,6406 


44,0807 


7,5984 


0,0422 


207 


5,248 


6,973 


6,7E7 


1,0011 E7 


1,0E8 


7,6418 


45,0045 


7,6188 


0,0230 


208 


5,248 


6,973 


4,0 E7 


9,0848 E7 


9,1 E7 


7,6415 


45,0045 


7,6168 


0,0247 


209 


5,007 


6,980 


1,7E9 


3,1963 E9 


3,2 E9 


7,6444 


47,1704 


7,6274 


0,0170 


210 


5,439 


6,966 


3,6 E6 


1,2001 E7 


1,2E7 


7,6379 


43,4243 


7,6208 


0,0171 


212 


8,985 


7,380 


1,3 E-7 


3,0038 E-7 


3,0 E-7 


8,1658 


26,2887 


8,1227 


0,0431 


213 


8,567 


7,394 


2,4 E-6 


4, 1949 E-6 


4,2 E-6 


8,1706 


27,5711 


8,1403 


0,0303 


214 


7,865 


7,417 


1,1 E-4 


1,6014 E-4 


1,6 E-4 


8,1789 


30,0316 


8,1053 


0,0736 



Table 3: Radial coordinates of the most probable formation of the a-particle before a-decay of nuclei Po (with 
charge number Z — 84) at different mass numbers A (for each nucleus we use: i?min = 0, 11 fm, i?niax — 100 fm; 
in region from i?min to 5 fm number of intervals is 2000, in region from 5 fm to 7 fm number of intervals is 
7000, in region from 7 fm to i?max number of intervals is 1000; Ar = R2 — Rtom', values Q, i?BMP, Tcxp are 
used from [13] ; twkb is half- live determined by approach WKB in [13]) 

nucleus is essentially closer to the experimental Tcxp in a comparison with half-live twkb obtained in [13] by 
approach WKB. So, method MIR is stable in determination of the coordinate of formation of the a-particle and 
half-live. For Ar we obtain: 

0,01fm< Ar < 0,07fm. (93) 

Also we have analyzed possible values of the coordinate of formation of the a-particle for these nuclei 
concerning the potential in the form DI. It turns out that at different locations of the starting point the 
calculated half-lives by approach WKB are some farther to their experimental values. 

6 Conclusions and perspectives 

In this paper the method of multiple internal reflections in description of the a-decay of nucleus in the spherically 
symmetric approximation is presented. In the approach MIR the formalism of calculation of the amplitudes of 
wave function, described moving of the a-particle from the internal region outside with its tunneling through a 
realistic radial barrier of arbitrary shape, has been constructed at first time. We establish (at first time): 

• The method MIR gives convergent values for the amplitudes of transmission and reflection, the coefficients 
of penetrability and reflection, obtained in description of leaving of the a-particle relatively a potential with 
barrier in form of a number of rectangular steps (multi-steps potential), with tending of this potential 
to the realistic a-nucleus potential. Therefore, limit values of the amplitudes and coefficients can be 
considered as the exact values concerning the realistic a-nucleus potential (without application of WKB 
approximation) . 

• Error in determination of the coefficients and amplitudes by the method MIR concerning the realistic 
a-nucleus potential is determined by accuracy of description of this potential by the multi-steps potential 
and accuracy of computer calculations (in all calculations, we have |T + i? — 1| < 1, 5 • 10^^^). 

In approach MIR a dependence of the penetrability coefficient and half-live on the location of the starting point 
^^stait, from where a-particle begins to move outside, is observed. According to this property: 

• We define the radial coordinate of the most probable formation of the a-particle inside nucleus before its 
a-decay as such starting point, at moving of the a-particle from which outside the half-live, calculated by 
approach MIR, is maximally closed to its experimental value. 
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• For a number of nuclei with the same charge number Z = 84 we have obtained essentially closer values 
for half-live, calculated by such approach, to their experimental values in a comparison with half-lives 
obtained in fT3] by WKB approach. This result is stable for all studied nuclei. 

• Inaccuracy of determination of half-lives for arbitrary studied nucleus on the basis of WKB method 
in a comparison with its experimental value can be explained (first of all) by reduction of the above- 
barrier external part of the radial potential starting from the external turning point, while inclusion of 
such contribution (by the method MIR, up to 50-60 fm) improves essentially the agreement with the 
experimental half-live. 

One can consider the pointed above dependence of half-live on the starting point by approach MIR as a new 
method of determination of the radial coordinate of the most probable formation of the a-particle before the 
a-decay of the nucleus on the basis of analysis of the experimental half- live for this nucleus. 

After construction of the algorithm of exact determination of the amplitudes of wave function relatively 
the multi-steps barrier (where these steps have arbitrary shapes), the method of multiple internal reflection 
transforms into a power tool for determination of stationary and (non-stationary) characteristics with high 
accuracy (and without application of WKB approximation) in variety of spherically symmetric problems of 
scattering and decay. Being the method with exact solutions, it allows to study deeper unusual and amazing 
properties of quantum systems and potentials, that points out availability of further improvements of the 
method, and such researches can have physical sense. So, it can be interesting by use of this method to 
study a reflectionless property, which is most evidentally shown in simplest soliton-like potentials [22] (see also 
[531 [211 US] ) and shape invariant potentials [551 [57], to clarify whether the effect of reflectionless tunneling is 
possible in real (hermitian) radial potentials |28l [29]. We suppose that this method should open a way for 
construction of new types of space asymmetric deformations of ID reflectionless potentials |30j . will give new 
variants of angular deformations of 3D reflectionless potentials obtained in frameworks of group theory [31] [32] . 
A generalization of the method MIR into complex (non-hermitian) potentials should allow to study properties 
of barriers in scattering problems with these potentials studied by approaches of non-linear supersymmetry 
(for example, see |3l [Ml [SB [36l [37] ) and N-fold supersymmetry [Ml [3S [40l [HI [42] , by approaches of hidden 
supersymmetry [43l[44]. In application to such problems, the method MIR can be a good test of early obtained 
results on the basis of methods of SUSY QM, group theory and inverse problem approach. 
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